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Cross Reference To Related Application 

Reference is made to and priority claimed from U.S. 
provisional application Ser. No. 60/251,191, filed Dec. 4, 2000, 
entitled METHOD AND APPARATUS FOR DETERMINING COMPONENT FLOW 
RATES For A MULTIPHASE FLOW. 

Field of the Invention 

The present invention concerns techniques for the 
measurement of the flow rate of each component of a multiphase 
fluid. More particularly, the present invention concerns 
techniques based on measurements at different locations along a 
conduit for determining the flow rates of each component of a 
multiphase fluid within the conduit. 

Background of the Invention 

In extracting oil and gas from a formation, it is 
advantageous to monitor the flow rates of the different 
components of the production fluid, usually gas, oil and water. 
Such a three -component mixture presents a problem that is often 
more difficult to solve than that of determining the flow rates 
for binary mixtures. For a broad class of well mixed binary 
mixtures flowing in a conduit of known cross sectional area, for 
which the density and the sound speed of each component are 
known, measuring the sound speed of the mixture and the 
convection velocity of the mixture provides sufficient 
information to determine the flow rates of each of the 
components. The sound speed of the mixture is directly related 
to the phase fractions of each component in the mixture. More 



specifically, to determine the four unknowns, the phase fraction 
of each component and the flow rate of each component , four 
equations are used. One equation expresses that the sum of the 
component (phase) fractions is unity. Another is that the speed 
5 of sound in the mixture depends on the speed of sound and 

density of the individual components and the phase fraction of 
each constituent. Finally, there are two equations relating 
the flow velocities of the individual components to the overall 
mixture flow rate. 

10 The approach of using a single sound speed measurement and 

p a mixture flow rate measurement used in determining the 
~ component velocities in the case of a binary mixture, however, 
■M cannot be extended to solve the problem of finding the component 
H.* velocities of a fluid with more than two components. Additional 

15fy information is required beyond what is provided by a measurement 
of the speed of sound in the mixture and a measurement of the 
M= flow rate of the mixture. 

nj To provide the required additional information, the prior 

g art teaches multiphase flow meters that typically rely on 

2 0M= several so called orthogonal measurement systems (said to be a 

orthogonal because each measurement system provides information 
that is at least partially independent of the information 
provided by the other measurement systems) . Multiphase flow 
meters according to such an orthogonal system approach include 
25 meters based on multiple -energy nuclear sources, ultraviolet 

measurements, capacitance measurements, venturi effect 
measurements, and infrared measurements. 

An alternative approach to determining the component flow 
velocities of a multiphase fluid in a conduit is to determine 

3 0 the additional required information from multiple point 

measurements, i.e. from measurements of the same information at 
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different places along the conduit. For example, the speed of 
sound in the fluid, the flow rate of the fluid, and the pressure 
and temperature of the fluid would be made at two or more 
locations along the conduit. In addition, a multiphase flow 
5 model is used to relate the values of the parameters that are 

measured at one location to those measured at another location. 
These relationships can provide the additional constraints 
required to solve for additional parameters. Typically, the 
equations are nonlinear in several variables. Many methods can 
10 be employed to "solve" for the flow parameters of interest. 

One class of methods defines an error function based on the 
-calculated values of the parameters compared to the measured 
values of the parameters. The flow-related parameters sought 

j=Jj (i.e. the component velocities or component phase fractions) are 

SI 

15ry adjusted -iteratively until the error is minimized. In this 

yf context, the flow parameters that result in minimizing the value 

y, of error function are assumed to be correct. 

u 

t 

One company, previously Loke of Norway (now owned by 



FMC/KOS) , is considered by many to have pioneered the general 



DO 
Q 

20U approach based on multiple point measurements for production 
allocation measurements in oil and gas production facilities. 
One implementation by Loke of the general approach is their 
software called Idun, which uses conventional pressure and 
temperature measurements along with a choke position measurement 

25 and knowledge of the fluid property characteristics to estimate 

the component flow rates (or phase fractions) . The overall 
accuracy and robustness of the Idun approach is directly 
influenced by the type and quality of sensors available. 

Another approach to determining component flow rates in a 
3 0 multiphase fluid based on making measurements at multiple 

locations along a conduit carrying the multiphase fluid is that 
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based on a gradio-venturi system, which includes a venturi meter 
and employs a remote pressure sensor located several hundred 
feet above the venturi. The pressure difference between the 
pressure at the venturi and that at the remote transducer can be 
5 related to the flow rate and composition through a multiphase 

flow model, and can be used in conjunction with the pressure 
difference due to the flow thorough the venturi to estimate the 
component flow rates . Such an approach has several drawbacks . 
It requires a venturi, which is intrusive, and it has an 
10 accuracy limited in two-phase flow to ±10% of the total flow. 

Moreover, the accuracy degrades substantially in the presence of 
any significant entrained gas. 



What is needed is a system of measuring the component 

velocities and phase fractions of a fluid that includes at least 

y = 

15fy three components, and that is nonintrusive, sufficiently 

accurate, and that does not provide spurious solutions because 

M= of insufficient information. 

f=* 

Summary Of The Invention 

ib si 

□ 

ill Accordingly, the present invention provides an apparatus 

20 and corresponding method for determining component flow rates of 

a multiphase fluid in a conduit, the fluid consisting of at 
least three known components, the method including the steps of: 
measuring at each of two different positions along the conduit 
at least four mixture quantities; providing a speed of sound in 
25 each of the components at the measured pressures and 

temperatures; providing a trial value for each of either the 
component flow rates or the phase fractions; using a 
predetermined model to calculate values for the measured mixture 
quantities based on the trial values for each of either the 
30 component flow rates or the phase fractions; using a 
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predetermined error function to determine an error value; and 
using a predetermined optimizing algorithm to determine whether 
the calculated values are acceptable, and, if they are not, to 
provide a new trial value for each of either the component flow 
5 rates or the phase fractions. 

In a further aspect of the invention, the error function is 
the sum of the squares of the difference between the measured 
and calculated values at each point . 

In another, further aspect of the invention, the four 
10 mixture quantities are the sound speed, the flow velocity of the 

y multiphase fluid, the pressure and the temperature. 

yp 

i.n Brief Description of the Drawings 

■crsr 1. 1 . ...I. — -■ i I . . 

m 

ffi The above and other objects, features and advantages of the 

k= invention will become apparent from a consideration of the 

15 3 subsequent detailed description presented in connection with 

y, accompanying drawings, in which: 

fg Fig. 1 is a schematic/ block diagram of a system for 

<=J measuring component flow velocities for a multiphase fluid in 

a conduit, according to the present invention; 

20 Fig. 2 is a set of plots showing the results of a 

simulation of use of the present invention in case of no free 
gas being present at an upstream position on a conduit; 

Fig. 3 is a set of plots showing the results of a 
simulation of use of the present invention in case of 20% (a 
25 phase fraction) free gas being present at an upstream 

position on a conduit; and 

Fig. 4 is a set of plots showing the results of a 
simulation of use of the present invention in case of 40% 
free gas being present at an upstream position on a conduit. 
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Best Mode For Carrying Out The Invention 



The present invention for determining component flow rates 
(or, equivalently, phase fractions) for a multiphase fluid in a 
conduit is based on making measurements at different points 
5 along the conduit, as indicated in Fig. 1. The term flow rate 

is used here in reference to a point along a conduit to mean 
mass per unit time flowing past the point along the conduit. 
The flow rate at a point is related to the average flow velocity 
at the point by the cross sectional area of the conduit at the 
10 point. The term phase fraction of a component of a multiphase 

O fluid in a conduit (at a point along the conduit) is used here 
jjri to indicate the fraction of the total mass of a sample of the 

ijj fluid (at the point along the conduit) that is due to the 

y = 

jjj presence of the component . Knowing the flow rate of each 

15 f Jf component is equivalent to knowing the phase fractions and the 

* mixture flow rate. The flow rate for a component, at a point 

f7 along a conduit, is simply the product of the component flow 

fy velocity, the phase fraction for the component, and the total 

:S mass per unit length at the point along the conduit. 

fcaf* 

2 0 s In the particular embodiment of a system according to the 

invention indicated in Fig. 1, four mixture quantities, namely 
the sound speed, the (mixture) flow velocity, the pressure and 
the temperature, are measured at two different positions along 
the conduit 11, at position A and at position B, using a 

2 5 measurement module 12 at A and a measurement module 13 at B, 

respectively. In the preferred embodiment, it is assumed that 
the flow can be modeled as time stationary so that measurements 
at the two positions at the same time are considered to be 
related. In other embodiments, however, the assumption that the 

3 0 flow is time stationary can be relaxed. In such embodiments, 

measurements at position B (downstream of position A) would be 
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assumed to be related to measurements at position A only after a 
time interval depending on the flow velocity measured at 
position A, and possibly also at position B. 

Still referring to Fig. 1, according to the invention 
5 values of the measured quantities are provided to a multiphase 

flow modeler 15, which also receives _as input the speed of 
sound in each of the individual components. Such component 
speed of sound values are determined based on pressure and 
temperature measurements at positions A and B and on an assumed 
10 PVT model. The multiphase flow modeler 15 then relates the 

% measurements at position B to the measurements at position A, 
yy according to the methodology described below, which involves 
jSf using successive trial values of the component flow rates (or 

01 the phase fractions) and an error function that compares the 

fli 

15ffj measured mixture quantities with calculated mixture quantities 

2 based on the trial values, and iterating until the value of the 
L error function is acceptably small. In the preferred embodiment 

E 

ssji : 

of the invention, the system includes an error function 

ffl 

□ evaluator 16 that computes, for each new trial value of the 
20^ component quantities (component flow rates or phase fractions) , 
each new corresponding value of the error function. An 
optimizer 17 then determines whether the calculated values are 
close enough to the measured values based on a predetermined 
tolerance, and if not, determines new trial values for the 
25 component quantities. The component rates (or the phase 

fractions) that minimize the error are ultimately provided by 
the optimizer 17 as the output of the system. 

For an N-component mixture in a conduit, the component 
densities p t , phase fractions </> i , and component sound speeds c ( . 
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are related to the mixture sound speed c mix and mixture density 
Pmix by tne equations, 



N 

Pmix = Z tlPi > 



i = l 



1 f 4, 

P C 1 TTf £> C' 

r mix mix l ~ 1 rii 

N 

I>, = 1 

/=1 



7, and (2) 



(3) 



For fluids in a flexible conduit, the propagation velocity 
J of sound, or any other compressional wave, is influenced by the 
0"i structural properties of the conduit. For a fluid in a conduit 
fri that is a vacuum-backed tube of radius R , thickness t and 
10[_ = . Young's modulus E , the measured mixture sound speed c meflJ is 

jTj related to the infinite domain sound speed c mix and to the 
y= structural properties of the tube via the relation, 

^ 1 1 2R 

pc pc Et 

r mix meas r mix mix 

However, as mentioned above, for a mixture of three or more 
15 components (called here a mult i -component mixture, as opposed to 

a binary component mixture) , the speed of sound no longer 
uniquely identifies the phase fractions of the different 
components. Additional information is required. 

The above relations combine to provide a relation between 
2 0 measured sound speed c meas , input fluid and structural properties 

_2R 

represented by <J ~~TT < and unknown phase fractions q>. (where 

Et 
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the index refers to one or another of the components) , a 
relation that can be expressed in the form, 



= Pn 



P C 

rid. 



(5) 



where the index a refers to either position A or position B, 
and the index i refers to the different components of the fluid. 



if! 



151 



20 



25 



Specializing the invention to a production fluid containing oil, 
gas, and water 

The invention will now be described in an embodiment that 
is particularly advantageous for determining the component flow 
rates (or, equivalently , the phase fractions) for a production 
fluid being pumped from a formation containing oil, gas, and 
water. It should however be understood that the invention has 
more general applicability; as will be clear, nothing about the 
invention restricts it to fluids of oil, gas, and water. 

In the particular case of a three-phase (three -component) 
fluid consisting of gas, oil, and water, equation (5) becomes, 



1 



= Pn 



mix, a 



go 



2 

ga~ga 



+ 



TO 



p„„c 



P c 

r oa oa 



+ 



2 

wa"wa 



(6) 



with figa + (j> oa +$ wa ~ 1 , and where the index a refers to either 

position A or to position B along a conduit bearing the 
production fluid. 

The invention imposes, as one condition, mass flow 
continuity of the individual phases, and uses a model for the 
evolution of gas from the oil. As the production fluid is 
pumped from position A up to position B, closer to the surface, 
gas escapes out of the oil/gas solution. The present invention 
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uses a model that calculates the amount of gas at position B 
based on, among other factors, the amount of gas and oil at 
position A. Conceptually, the model tracks the multiphase flow 
as a four-component mixture, a mixture consisting of water, oil, 
5 gas present at position A, and gas evolved from the oil between 

position A and position B. 

Application to well-mixed flow 

■ The invention will now be described in the particular case 
of well -mixed flow of a fluid consisting of oil, gas, and water. 
10^ The extension to the general case, i.e. to the case of flow that 
yg is not well-mixed, will be apparent to one skilled in the art. 

yg In the case of well -mixed flow, there is no slip among the 

'^2 individual phases, i.e. the mixture is homogeneous. In such a 

fU mixture, the condition of mass flow continuity from position A 

15 U5 to position B becomes, 

H- Pmix,A A A U A = Pmix,B A B U B , (7) 

JS[ where A 0 is the cross section of the conduit at the position 

indicated by the subscript a , and u a is the corresponding flow 

rate. In terms. of the individual phases, assuming that the 
20 cross sections at position A and position B are equal, the 

condition of mass flow continuity expressed by equation (7) 
becomes , 

U A (toAPoA + tgtPgA + <f> W APwA ) = U B [toBPoB + ^PgB + </>„ B Py»B ) • ( 8 ) 

Conservation of the individual phases results in three 
25 equations, one equation for each of the three components. For 

the oil and water components, the invention in the preferred 
embodiment assumes that the density of the water and oil 
components is the same at position A and position B, yielding, 
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10% expansion characterized by its compressibility factor. 

%£J Secondly, it is assumed that amount of gas that evolves from the 

fluid as it moves from position A to position B is proportional 
r=J to the amount of oil present at position A. Expressing these 
relations in terms of volumetric flow rates at standard 
15M= temperature and pressure (indicated by the subscript R)_ yields 
l=li the following relations, 



U A^oAPoA = U B</>obPoB, (9) 

and 

u a0 w aP w a = u b</> w bP w b ■ do) 

These simplifying assumptions can be relaxed without altering 
the basis of the approach of the present invention. 

The third equation derives from the conservation of the gas 
phase. The present invention models two mechanisms to account 
for the change in volumetric flow rate of the gas. First, it is 
assumed that the gas present at position A under goes an 



y u A ^A A ^^^- + fi R (P B -P A )(/> oA u A A A =uJ gB A B ^-^-^- (ii) 

^ A "r *A *R *B 

where /3 is the so-called gas evolution constant, which is 
defined as the amount of gas (in standard cubic feet) evolved 
2 0 from a barrel of oil as the pressure is lowered one pound per 

square inch (psi) , and Z is the gas compressibility factor, 
defined as, 

VPT 

Z^. (12) 
V R P R T 

The gas evolution constant J3 can be determined from a standard 
25 differential vaporization test in which the amount of dissolved 
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gas (in standard cubic feet) per barrel of residual oil is 
reported as a function of pressure (see e.g. PVT and Phase 
Behavior of Petroleum Reservoir Fluids, by Ali Danesh, Elsevier, 
1998) . The gas compressibility factor Z is also reported in 
typical PVT field analyses of black oils. 

The first term in equation (11) represents the amount of 
gas present at position A, and the second term represents the 
amount of gas evolved from the oil as the fluid moves from 
position A to position B. 

The invention therefore constructs a system of seven 
equations. The sound speed mixing laws, given by equation (6), 
provide two equations. The phase conservation conditions, given 
by equations (9-11), provide three equations. These seven 
independent equations are in terms of six unknowns (the 
component flow rates or the phase fractions at both positions A 
and B) and so provide more than enough information to uniquely 
determine the six unknowns. Since the system of equations is 
over determined, instead of attempting to solve the system of 
equations (which may include inconsistent conditions) , the 
invention seeks a solution that minimizes an error function. In 
the preferred embodiment, the invention uses as an error 
function a simple sum of the squares of the differences between 
the calculated values of a quantity measured at the two 
positions A and B and the measured values, i.e. using F i to 
represent a quantity that is measured at position A and position- 
B, where the index i ranges over N such quantities, the 
invention uses as the error function, 

N i 2 

E ~ ^ Ycalcj ~ Fmeasj \ (13) 



The application of this general approach can take many forms, 
i.e. many different algorithms or procedures may be used to 
minimize the error function E ; how, specifically, the error 
function is minimized is not the subject of the present 
invention. 

The accuracy of the general approach used in the present 
invention, that of minimizing an error function as explained 
above, depends on many factors, including: the fidelity of the 
parametric model; the quality, quantity and types of 
measurements; the accuracy of the calibration of the model input 
parameters; the accuracy of the well parameters, including for 
example the radius R , thickness t and Young's modulus E of the 
conduit; and the accuracy of the fluid parameters, including for 
example the component speeds of sound. 

The inventors conducted several Monte Carlo simulations to 
evaluate the multiphase flow measurement approach developed 
above. The well parameters were selected to be representative 
of a producing black oil. The parameters used in the simulation 
are listed below. 

Table 1. Parameters of components. 



Fluid 


Density 
(kg/m"3) 


sound speed 
(ft/sec) 


oil 


700 


3500 


water 


1000 


5000 


gas 


1.2*P/P atra 


1200 



Table 2. Parameters of the mixture and conduit. 



Component- independent parameters 


Pressure at A 


3000 psi 


Pressure at B 


2800 psi 


Temperature at A and B 


100 C 


p, gas evolution constant 


0.316 scf/stb/psi 



Inner diameter of conduit 


3 . 0 inches 


Wall thickness at A 


0.50 inches 


Wall thickness at B 


0.25 inches 



Figs. 2-4 present the results of the Monte Carlo 
simulations, in which the input parameters and measured 
parameters are varied in a random, but bounded, manner. The 
three figures illustrate the ability of the invention, even 
5 using the simplifying assumptions noted above, to measure the 

multiphase flow parameters as a function of oil and water cut 
for three levels of free gas present at position A (namely 0% 
O for Fig. 2, 2 0% for Fig. 3, and 4 0% for Fig. 4) in the presence 
.}=! of the uncertainty for each of the parameters indicated in Table 

lOijO 3 below. The pressure was not directly varied in the 

oi 

fj\ simulations; the effect of any uncertainty in pressure is 

^ accounted for in the variations of the other model parameters. 

f Table 3 . Uncertainty in model parameters . 



Parameter 


Uncertainty 


P, gas evolution constant 


2% 


Oil sound speed 


2% 


Gas sound speed 


1% 


Water sound speed 


1% 


Oil density 


1% 


Gas density 


2% 


Water density 


0.5% 


Measured sound speed at A and B 


0.25% 


Measured flow velocity at A and B 


0.25% 



For each of the simulations indicated in Figs. 2-4, a single 
15 wall sub (measurement station) is assumed, with a thickness at 

position A of 0.5", and a thickness at position B of 0.25", as 
indicated in Table 2 . 
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Discussion 

The simulations indicate that the invention provides 
reasonable estimates of oil, gas, and water phase fractions, 
within ±10% of reference values for a reasonable level of 
uncertainty in both measured and input parameters for free gas 
volume fractions ranging from 0% at position A to 40% at 
position A. 

The solution appears to be robust to uncertainty; gradually 
increasing the uncertainty in the input and measured parameters 
10_, appears to translate into only a gradual degrading in flow 
measurement accuracy. 



yl Scope of the Invention 

On 

py It is to be understood that the above -described 

m 

arrangements are only illustrative of the application of the 
15M-- principles of the present invention. Numerous other 

modifications and alternative arrangements may be devised by 
those skilled in the art without departing from the spirit and 
scope of the present invention, and the appended claims are 
intended to cover such modifications and arrangements. 
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